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Abstract 

We show that the Auger Air Shower Array has the potential to detect neutri- 
nos of energies in the 10 19 eV range through horizontal air showers. Assuming 
some simple conservative trigger requirements we obtain the acceptance for 
horizontal air showers as induced by high energy neutrinos by two alternative 
methods and we then give the expected event rates for a variety of neutrino 
fluxes as predicted in different models which are used for reference. 

PACS numbers: 95.85.Ry, 96.40.Tv, 96.40.Pq, 98.70.Sa 
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I. INTRODUCTION 



High energy neutrino detection is one of the experimental challenges in particle astro- 
physics for the forthcoming years because it opens a new window to the regions of the Uni- 
verse that are otherwise shielded from us by large amounts of matter. It is widely believed 
that one of the most appropriate techniques for neutrino detection consists of detecting the 
Cerenkov light from muons or showers produced by the neutrino interactions in underground 
water or ice. This allows the instrumentation of large enough volumes to compensate for 
both the low neutrino cross section and the low fluxes expected. There are several projects 
under way to build sufficiently large detectors to measure the expected signals from a variety 
of sources Q . 

On the other hand many years ago it was suggested that deeply penetrating high energy 
particles, such as muons and neutrinos, initiate large horizontal air showers that can be 
detected at ground level || . At large zenith angles the electromagnetic part of ordinary air 
showers, initiated by cosmic rays of hadronic (or electromagnetic) nature, is attenuated by 
the atmosphere well before reaching ground level. This should allow the identification of the 
showers initiated by these deeply penetrating particles. 

Recently there has been a proposal to build two 3000 km 2 extensive air shower arrays, 
one in each hemisphere, to detect cosmic rays with energy above 10 19 eV (the Pierre Auger 
project) [||. Each array consists of a hexagonal grid of water tanks, 3.5 m diameter and 1.2 
m height, combined with an optical fluorescence detector. The tanks are separated 1.5 km 
from each other and instrumented with photodetectors to detect the Cerenkov light emitted 
as photons and charged particles from the shower front cross it. We will show that these 
detectors can also play an important and complementary role for detecting neutrinos of 
energy around 10 19 eV and above. 

Neutrinos of these energies cannot penetrate the earth so any detector searching them 
must look for events with zenith angles between near horizontal and vertical downgoing. 
For showers of sufficient energy the array efficiency is high and the low target density for 
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neutrinos provided by the atmosphere is compensated by the huge surface area of the planned 
array. The observatory will complement other neutrino detectors in construction or planning. 

In this article we show that the Pierre Auger Observatories have an enormous potential 
to detect neutrinos in the EeV range. In section II we firstly present some of the possible 
EeV neutrino fluxes that have been discussed in the literature which will be used later in 
the evaluation of the event rates expected by the Auger particle arrays. In section III we 
discuss the acceptance for the horizontal air showers. Firstly we obtain a rough but intuitive 
analytical expression for the acceptance of a large array, which we discuss in the context of 
the Pierre Auger project. In section III. A we present an acceptance calculation for the Auger 
Project based on a geometrical approach. In section III.B an alternative and conservative 
calculation based on Monte Carlo simulation is presented. In section IV we estimate the 
horizontal shower rates expected for the neutrino fluxes addressed in section II, and section 
V is reserved for the conclusions. 

We do not attempt to discuss event reconstruction or backgrounds which will be ad- 
dressed elsewhere. 

II. EEV NEUTRINO FLUXES 

Neutrinos of EeV energies are likely to be directly produced together with cosmic rays 
of the highest energies whose origin is still a matter of speculation. Several possible sources 
of neutrinos in the EeV region have been discussed in this context. Moreover, provided 
that the highest energy cosmic rays are extragalactic in origin, as it is currently believed, 
EeV neutrinos have to be produced in their interactions with cosmic microwave photons. 
We will use some of the neutrino flux predictions reported in the literature as reference 
calculations to evaluate the event rates that could be expected in each Auger Observatory. 
In this section we briefly motivate these fluxes which are shown in Fig. 1 compared with the 
expected atmospheric neutrino flux. 

Active Galactic Nuclei (AGN), the most energetic objects known, emit most of their 
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luminosity in gamma rays J3J and may be the source of the highest energy cosmic rays 0. 
Multiwavelength observation of these sources leads us to believe they are powered by large 
accreting black holes, where shock fronts provide particle acceleration. If, besides electrons, 
protons are also accelerated, as some models suggest, photoproduction of pions with the 
ambient light plays a crucial role (see for example Ref. || and references therein). The 
energetic gamma rays observed come from the decays of 7r°'s and high energy neutrinos 
from the decays of charged pions become a signature of these models. First models for 
neutrino production in AGN assumed shock acceleration in the AGN cores and predicted 
relatively flat fluxes up to energies of about 10 15 eV. For our event rate calculation we 
select the prediction of Ref. [0], labelled AGN-92C, which is quite similar to that in Ref. 
(AGN-95C). 

There is however recent evidence that the GeV to TeV gamma ray emission observed 
from AGN corresponds to the blazar class 0. In a unified AGN description, the blazar class 
is identified as AGN with jets of ultrarelativistic particles streaming out of their cores and 
pointing towards us. Most recent models for the proton blazars site the acceleration in the 
jets themselves. Photoproduced neutrinos are Lorentz boosted to energies well above those 
predicted in the models of acceleration in the AGN cores. We use the prediction of Ref. || 
(labelled AGN-95J) which illustrates that the emitted neutrinos may extend well into the 
EeV region in agreement with Ref. ||10|| . 

We also consider the more speculative and uncertain models where the highest energy 
cosmic rays are produced by the decays of topological defects. These objects emit massive 
X particles, predicted in Grand Unified theories, which decay and fragment into Standard 
Model particles. There are a fair number of neutrino flux predictions from these models 
depending on the different time evolution of the effective injection rate of X particles per 
unit volume (usually denoted with a parameter p) , the masses of the X particles themselves 
and other unknown parameters. Some models have already been constrained by bounds 



on horizontal showers induced by neutrinos |TT[ and also by measurements of the 100 MeV 



diffuse gamma ray background |L2[ . We select two fluxes illustrating the range of predictions 



in such models. We use the prediction with p — 1.5 in Ref. |L3| (labelled TD-92) and that 



given in Ref. [|14[] which is labelled as TD-96. TD-92 is the lowest neutrino prediction given 
in Ref. |13|] and is not severely affected by experimental constraints. 

Lastly we consider two predictions for the high energy neutrinos produced in the inter- 
actions of cosmic rays with the cosmic microwave background: CMB-91 [|TJJ and CMB-93 
corresponding to model 4 in Ref. | 16|| . They illustrate the variations that can be associated 
to various uncertainties intrinsic to these calculations. These neutrinos are a direct result 
of the Greisen-Zatsepin-Kuz'min cutoff and must be produced if the highest energy cosmic 
rays are of extragalactic origin. Unfortunately it is possible that the flux levels are too low 
to be observed either with the Pierre Auger Observatory or with other neutrino detectors 
in planning or construction. 

When the electron neutrino flux prediction is not explicitly available in any of the models 
used, we approximate it to be a factor of two below the muon neutrino flux as expected by 
naive channel counting in pion production and decay. 



III. ACCEPTANCE FOR NEUTRINO SHOWERS 

Neutrinos produce showers in most interactions with the atmosphere which are of dif- 
ferent nature depending on the process in consideration. We consider both deep inelastic 
charged and neutral current interactions which always produce hadronic type showers. In 
the case of charged current electron neutrino interactions the emerging electron contributes 
in addition a pure electromagnetic shower carrying a large fraction of the incoming particle 
energy. We will ignore the resonant cross section because it is only significant near the peak 
of the cross section which occurs at an incoming neutrino energy of 6.4 PeV, well below the 
region of high efficiency for the Pierre Auger project. 

For a neutrino flux d§ v /dE v interacting through a process with differential cross section 
da/dy, where y is the fraction of the incident particle energy transferred to the target, the 
event rate for horizontal showers can be obtained by a simple convolution: 
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$ sh [E sh > E th ] = N aPair / dE sh / dy -r^{E v ) —(E u , y) A(y, E v ) (1) 

J E th Jo dh/ v ay 

where N a is Avogadro's number and p air is the air density. The energy integral corresponds 
to the shower energy E s h which is related to the primary neutrino energy E v in a different 
way depending on the interaction being considered. A is a geometric acceptance, a function 
of shower energy, which corresponds to the volume and solid angle integrals for different 
shower positions and orientations with respect to the array. The function is different for 
showers induced by charged current electron neutrino interactions from those arising in 
neutral current or muon neutrino interactions. This is because hadronic and electromagnetic 
showers have differences in the particle distributions functions, particularly for muons. 

For showers of sufficient energy, incident with zenith angle 9 zenit h, the effective volume 
for horizontal shower detection is given by S cos 9 zenit h, where S is the surface area covered 
by the array, multiplied by the range of allowed positions for the first interaction point along 
the incident direction. We can estimate this range approximating shower maximum as a 
disk of some effective radius r. This radius should approximately be the maximum distance 
to the shower axis at which the particle density is high enough to trigger the detector 
tanks (see below). As the first interaction point is moved along the incident direction the 
intersection of this disk with the detector plane is an ellipse with major axis q = 2r/sina 
with a = 90° — 6 zenith . The projection of this axis onto shower direction (q cosa) is the 
wanted range (see Fig. 2). 

We estimate A integrating this volume over the possible solid angle orientations of the 
shower, dQ = 2nd(sina), and restrict the integration to horizontal showers i.e. 0° < a < 
a max — 15°. Approximating cosa ~ 1 we obtain a simple analytic expression in terms of r. 

/ d(sma) sin a h / d(sma) sin a W) 

Jsinai sin a Jo / 

= S x 2n r (2 sin a max — sin cti) (2) 
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For a < ot\ = sin _1 (2r /W) the ellipse axis q exceeds W, the average size of the array Q. The 
last term in the brackets of Eq. (2) represents a small correction obtained when restricting 
the range of the interaction point to be below the longitudinal size of the array. Most 
importantly the acceptance is seen to depend on the array surface area and the shower 
radius. The estimate makes the assumption that showers with shower maximum intercepting 
the array are detected with close to 100% efficiency The effective radius has to be chosen 
to match this requirement which will only hold for energies above a given threshold. It 
is conservative in the sense that it ignores the fact that showers can trigger the detector 
without having shower maximum intercepting the array, or even without the shower axis 
going through the array (i.e. completely horizontal showers). 

It is easy to see that the acceptance for the Auger particle arrays is comparable to other 



neutrino detectors in planning p7| . Since design requirements have led to a tank size that 
allows near 100% efficiency for vertical showers of energy above 10 19 eV ||, the expected 
signal in the tanks should be large enough to trigger up to distances from shower axis of order 
the separation between the tanks, 1.5 km. It is thus conservative to expect similar efficiencies 
for horizontal showers above 10 19 eV as they should not differ that much from vertical ones at 
detector level and to approximate r ~ 1.5 km. The water Cerenkov technique is particularly 
well suited for horizontal showers. The transverse separation between detectors will be 
substantially reduced for near horizontal showers and the extra depth of the tanks in the 
horizontal mode should enhance the signal from the muons in the shower, compensating the 
reduction in transverse area of the tank to the incident direction. 

Taking r = 1.5 km {on = 4°) we obtain an estimate for the Auger array acceptance of 
showers of energy above ~ 10 19 eV, A = 13000 km 3 sr which when multiplied by an air 
density p a i r ~ 10~ 3 g cm~ 3 gives 1.3 10 7 kT sr. Underground high energy neutrino detectors 



*For S = 3000 km 2 the "diameter" of the Pierre Auger array is approximately D = 65 km and 
W ~ 0.70L> ~ 45 km. 
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in planning aim towards an active volume in the range of 1 km 3 of water [Tj|. In models 
where most of the events are due to neutrinos well above the PeV region, the Earth will be 
opaque to them. The corresponding acceptance for contained events is of order 6 10 6 kT sr[], 
illustrating how the Pierre Auger project may come into play. 

The above estimate of the acceptance does not display the energy dependence but it 
illustrates how it depends on the effective shower radius, which in turn must depend on 
energy. As it happens for vertical showers, the acceptance must show a kind of threshold 
behavior, since sufficiently low energy showers cannot be detected with such a sparse array. 
The dependence of the acceptance on shower energy has been introduced in two alternative 
ways. In a more geometrical approach we have followed on the above line of thinking consid- 
ering in more detail the shower depth development and the lateral distribution functions for 
electrons and positrons. As an alternative we have simulated showers at different positions 
with respect to the detector in a Monte Carlo approach. 



A. Geometrical Approach 

The calculation treats electromagnetic and hadronic cascades independently using the 
conventional NKG charged particle lateral distribution functions normalized with the 
Gaisser |20| (Greisen (T^j) parametrization for the total number of electrons and positrons 
in hadronic (electromagnetic) cascades. 

We assume that triggering requirements for each tank can be specified as fixed numbers 
for the electron density in the electromagnetic and hadronic showers: The relevant quan- 
tity for the tank signal is the Cerenkov light, proportional to the charged particle tracklength, 
which is in turn approximately proportional to energy deposition. Muons travel through the 
tank depositing about 2 MeV/(g cm -2 ) while most electrons and photons are stopped in 



^ However for muon neutrinos the long range of the muon increases the acceptance for non con- 
tained charged current interactions. 
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the tanks, depositing all their energy. We use p* because there are convenient parameteri- 
zations in the literature for this quantity. The relation between pf 1 and energy deposition is 
however dependent on the shower position and can only be done in an approximate manner. 

The above assumption naturally defines the active region of a shower, the volume within 
which the electron density exceeds the given threshold p* h . The region is bounded by a cigar 
shaped contour plot of the three dimensional electron distribution function at p 1 ^ 1 . For the 
solid angle- volume acceptance integration, we consider the number of tanks contained in the 
intersection of these active volumes with the detector plane as the shower directions and first 
interaction point are varied. The electron-positron density at these tanks is by definition 
above threshold. At typical large distances (of order a kilometer) to the shower axis, we 
estimate that an electron density in the range 0.6-1.2 m -2 (0.3-0.5 m -2 ) for horizontal 
electromagnetic (hadronic) showers is equivalent to an energy deposition of 500 MeV in a 
tank (corresponding to about two vertical muons). This is considered sufficient for shower 
detection ||. 

Approximating the intersections of these active regions as rectangles allows a complex 
but analytical solution to the problem if the only orientations and positions of the rectangles 
considered in the integration are those which contain at least n tanks in the same row (see 
Fig. 3). As a result we obtain the acceptance curves as a function of shower energy for 
each shower distribution function which only depend on the parameter The important 
advantage of this calculational method against the obvious alternative which involves Monte 
Carlo simulation (addressed below) is computing speed. The approach converts in straight- 
forward the otherwise lengthy evaluation of the effect of changing the input parameters 
in the calculation such as trigger conditions, tank thresholds, parameterizations of shower 
distribution functions and even array spacing. 

In Fig. 4 we show the acceptance results for a trigger of at least n = 3 aligned tanks, 
having the conservative electron density of 1.2 m -2 (0.5 m -2 ) for electromagnetic (hadronic) 
showers and considering only showers with zenith angle higher than 75 degrees. For each 
shower type there are two curves, one considers showers with axis falling within the array 
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and the second also includes an approximate calculation of the triggering showers with axis 
not going through the array. It is important to note that these showers make most of the 
contribution to the acceptance in the high energy limit (corresponding to huge showers flying 
parallel and on top of the array). 

Besides using pf- as the parameter for tank threshold, the geometrical algorithm de- 
scribed above has a potentially more important drawback. By using parameterizations for 
air showers, it ignores the effect of the ground on shower development. The simplification 
can be thought to grossly overestimate the acceptance. Using EGS4 with constant density 
air, we have simulated the effect of the ground for very inclined electromagnetic showers of 
energies up to 10 15 eV. In spite of the reduction in the lateral distribution observed after 
the core hits the ground, there is only a very small reduction in the acceptance assuming the 
results can be extrapolated to higher energies. In any case we have decided to contrast the 
acceptance results against a totally different calculational approach, which we chose to rem- 
edy these drawbacks in a conservative manner. That way we hope to bracket a prediction 
for the horizontal shower acceptance of the Auger project. 

B. Monte Carlo approach 

We have performed a totally independent and conservative calculation with a hybrid 
Monte Carlo and parametrization method. The approach uses the simulated curves of 
energy deposition versus distance from shower axis in the Pierre Auger water tanks for a 
10 19 eV shower at different depths. These results are then simply scaled with shower energy. 
We have restricted the range of depths for the shower to the interval [500,1600] g cm~ 2 and 
we only consider showers with the shower axis falling within the array. In our effort to be 
conservative we completely neglect the contributions from particles after the shower core 
hits the ground. We have then generated horizontal events by allowing the starting point 
and direction (75° < 6 zenit h < 90°) of the shower to vary in a random way, only restricted 
by the array size and a maximum height of 3 km. 
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In the simulation of the acceptance, events are detected when the shower deposits more 
than 500 MeV of energy in each of four or more tanks. We have used an average pa- 
rameterization assuming half of the energy induces a hadronic shower and the other half 
an electromagnetic shower. Besides correcting drawbacks in the geometrical approach, the 
method has the advantage of allowing the study of the different patterns of hit tanks on 



an event by event basis ||21|| . The results are shown as points in Fig. 4 for some discrete 
energies. Given the fundamental differences between the two approaches, it is not possible 
to test one against the other. Nevertheless it is fair to remark that the agreement is rather 
good and that the differences of results can be interpreted on the basis of the different in- 
puts intrinsic to each method. We thus take the Monte Carlo approach as a lower bound 
of the acceptance and the difference between the two approaches can be considered as an 
indication of the degree of uncertainty. 

For the detection of the horizontal showers included in the acceptance calculations de- 
scribed here, only directional reconstruction and an ability to separate the electromagnetic 
part of the shower from the muonic component should be sufficient in principle. Neverthe- 
less if the long range of the muons produced in the horizontal showers of hadronic nature 
is also taken into account, it can be argued that the acceptance may exceed that actually 
calculated in this work. The muon component effectively enhances the length of the active 
region of the shower, what must increase the acceptance for deeply penetrating particles. 
The identification of such showers would become however more complicated, because they 
look more like those induced by cosmic rays at very large zenith angles. 

IV. EVENT RATES 

To calculate event rates we use the deep inelastic charged and neutral current cross- 
section from two sets of structure functions MRS(Rl) and GRV (MS renormalization scheme) 



parton distributions [22]. With increasing energy, neutrinos interact with partons carrying 



a lower fraction x of nucleon momentum which extends beyond the kinematical limits of 
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the parton density parameterizations. For the first set we extrapolate to low x beyond 
validity of the parameterization, using the slope of xq(x), at the lowest x permitted (in this 
case x = 10~ 5 ), where q(x) is the standard parton distribution. This is consistent with 



extrapolations based on the leading log approximation [p^j . The second set, GRV, can be 
cautiously used on its own for low x. 

Table 1 displays the event rates for different neutrinos fluxes, calculated with Eq. (1) 
for the each of the two cross sections. The first entry corresponds to the result of our 
conservative approach based on the simulation and the higher one to the acceptance results 
calculated in the geometrical approach including showers with axis not intercepting the array. 
Recent models of proton acceleration in AGN jets (AGN-95J) and some models of decays of 
topological defects (such as TD-92) predict neutrino fluxes giving measurable rates in the 
most conservative assumptions. For neutrinos from AGN cores the majority of the detected 
showers lie in the 10 15 — 10 17 eV region. This corresponds to the threshold region for the 
acceptance curves where the efficiency is low and are very sensitive to trigger conditions. 
For these neutrino fluxes, a 1 km 3 conventional underground neutrino detector is expected 
to give more events. The detection of neutrinos from interactions of cosmic rays with the 
cosmic ray background is unfortunately only marginal for the most optimistic predictions. 
It should be mentioned that conventional muon underground detectors in planning will have 
similar difficulties, if not higher, in the detection of these neutrinos. 



V. CONCLUSIONS 

The Pierre Auger project can be made sensitive to ultra high energy neutrino fluxes 
through horizontal showers if an appropriate trigger is implemented. Its acceptance for 
detecting contained neutrino events of energy above E u ~ 10 19 eV can be very large and 
comparable to neutrino telescopes in planning. The highest efficiency for horizontal shower 
detection with the Pierre Auger Project is expected at energies about 10 19 eV, what makes 
it in principle a tool for the search of the neutrinos from interactions of cosmic rays with the 
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cosmic microwave background. The Auger array can detect the very high energy neutrinos 
from the decays of topological defects in some predictions. In any case the plethora of 
topological defect models will be further constrained by the neutrino detection capabilities 
of the Auger observatories. Most recent predictions for neutrinos produced by protons 
accelerated in the jets of AGN have also prospects to be detected as horizontal shower 
events in the Auger particle arrays. 

In any case the target mass of the Auger detector is sufficiently large to serve as a 
significant explorating tool for ultra high energy neutrinos, what ever their source. 
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TABLES 



v source 



AGN - 92C 
AGN - 95C © 
AGW - 95 J @ 
CM5 - 91 



CMB - 93 P] 



TL> - 92 P| 



- 96 114 



MRS(Rl) 



2.5/0.4 
1.5/0.2 

7/2.1 
0.5/0.2 
0.2/0.06 

12/4 
1.4/0.5 



GRV95 



2/0.3 
1.2/0.2 
5.5/1.6 
0.4/0.1 
0.13/0.04 

7/2.4 
0.9/0.3 



Energy (GeV) 



10 6 < E < 10 s 
10 6 < E < 10 s 
10 8 < E < 10 10 
10 10 < E < 10 11 

10 8 < E < 10 12 

10 9 < E < 10 13 
10 10 < E < 10 12 



TABLE I. Yearly neutrino event rates for different diffuse fluxes and cross sections (see text). 
The first entry corresponds to the upper set of curves in Fig. 2 (all showers) and the more conser- 
vative second entry is the result of the simulation. The last column gives the energy range of the 
bulk of the events. Events correspond to a single array of 3000 km 2 . 



Figure 1: Neutrino flux predictions in the EeV range as labeled in the text. 



Figure 2: Schematic picture of a horizontal shower, illustrating the disk of radius r at shower 
maximum and the cylinder it spans as the first interaction point is shifted along the incident 
direction. The intersection of the cylinder is an ellipse of major axis q (see text). 



Figure 3: Illustration of the intersection of the "active part" of the shower (as described in 
the text) and the detector plane used in the geometrical approach. The intersection, which 
is close to an ellipse, is approximated by a rectangle as illustrated. The dots correspond to 
the detector tanks. 



Figure 4: Acceptance of the Pierre Auger detector to near horizontal showers {6 zenit h > 75 
degrees). Volume units are km 3 of water equivalent. Crosses are for the results in the 
Monte Carlo approach, lines are correspond to the geometrical integration for electromag- 
netic showers (dashed) and hadronic (solid). The lower set of curves corresponds to only 
showers with axis falling in the array and the upper set takes all showers into consideration. 



15 




Fig. 3 



